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Abstract 

This paper describes the details of implemen- 
tation of a general numerical procedure developed 
for the accurate and economical computation of 
natural frequencies and associated modes of any 
elastic structure rotating along an arbitrary 
axis. A block version of the Lanczos algorithm is 
derived for the solution that fully exploits asso- 
ciated matrix sparsity and employs only real num- 
bers in all relevant computations. It is also 
capable of determining multiple roots and proves 
to be most efficient when compared to other, simi- 
lar, existing techniques. 

Introduction 

Gyroscopic structural systems are often 
encountered in practice. Thus, some structures 
such as satellites are usually spin-stabilized 
whereas others, like helicopters and turbines, 
have rotating parts. An accurate evaluation of 
their frequencies and mode shapes is of utmost 
importance in predicting their stability and also 
in implementing effective closed-loop control of 
the gyroscopic systems. The usual solution proc- 
ess starts with a finite-element discretization 
of the structure yielding appropriate stiffness 
and inertia properties, which in turn are uti- 
lized to yield the natural frequencies and asso- 
ciated modes. Such data are next utilized to 
compute unsteady aerodynamic forces enabling com- 
putation of flutter and divergence charac- 
teristics. An extension of the analysis yields 
the state-space matrices enabling open- and 
closed-loop control analysis of the structure. 

The accuracy of such an analysis is, however, 
entirely dependent on appropriate computation of 
vibrational characteristics of the system. 

The equation of free vibration of any struc- 
ture discretized by the finite-element method and 
spinning along an arbitrary axis with a uniform 
spin rate fl is given by 

WJ + Cq + Kq = 0 (1) 

in which 

K = K E + Kq + K 4 

and 

M inerti a matri x 
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C skew-symmetric Coriolis matrix, function 

of u 

K E elastic stiffness matrix 

Kq geometric stiffness matrix, function of fl 2 

K* centrifugal force matrix, also a func- 
tion of ft 2 

q deflection vector 

For small vibrations, the K and M matrices are 
real, symmetric, and positive definite. The 
solution of Eq. (1) may be achieved by first 
rearranging the same as 

Ay + By = 0 (2) 

in which 



where A is symmetric and B is skew-symmetric. A 
solution of Eq. (2) is obtained by substituting 
y = e wt yielding 

(A + u>B)y = 0 (3) 

in which the natural frequencies u> are pure imag- 
inary, the vectors being complex and both occur- 
ring as complex conjugate pairs. 

The conventional solution process for Eq. (3) 
involves implicit inversion of A to reduce the 
eigenvalue problem in terms of a single matrix 
of order twice the original size, which however 
is rather inefficient due to its nonsparse char- 
acter and increased order. 

A combined Sturm sequence and inverse itera- 
tion (SS/II) procedure was presented earlier 1 for 
the eigenproblem solution of gyroscopic systems 
that exploits inherent sparsity of constituent 
matrices of Eq. (2). Reference 2 provides a sur- 



vey of solution methods for free-vi bration analy- 
sis of structures including spinning ones. An 
improved version of the SS/II technique has further 
been presented in a recent paper 3 that also gives 
details of numerical techniques for computation of 
in- and out-of-plane forces in a shell and also 
line elements spinning along an arbitrary axis. 

The Lanczos method has been applied earlier 4 * 5 for 
the eigenproblem solution of real symmetric matri- 
ces, Reference 4 also presents the relative merit 
of the block Lanczos algorithm over the conven- 
tional nonblock procedure. In a recent paper, 5 a 
nonblock version of the Lanczos algorithm was pre- 
sented that is suitable for the economical solution 
of the eigenproblem of gyroscopic systems. 

The purpose of this paper is to provide details 
of a related block Lanczos algorithm and its imple- 
mentation in a general-purpose finite-element 
computer program, STARS 7 (STructural Analysis 
Routines). Numerical results are also presented 
that prove the efficacy of the current solution 
technique. 

Implementation of a Block Lanczos 
Solution Procedure 

To implement the current procedure, Eq. (3) 
is first rewritten as 

(A - AD)y = 0 (4) 

in which D = i*B is a pure imaginary Hermitian 

matrix, i* is the imaginary number V^T, the roots 
A = i*w are real and occur in pairs X\, -Aj, ..., 

A n , - A n whereas the corresponding eigenvectors 
occur in complex conjugate pairs. The roots of 
the original system defined by Eq. (3) may then 
be simply obtained as u = A/i* while noting that 
the eigenvectors remain the same for both cases. 

To develop the present algorithm, it is nec- 
essary to yield a single matrix out of the set 
of two matrices that define Eq. (4). This is 
achieved by performing a Choleski decomposition 


It may be noted that H is a pure imaginary 
Hermitian matrix and the current transformation 
retains the banded form of associate matrices. 

In all subsequent computations, n defines the 
order of H, whereas m^ denotes the half-bandwidth 
of constituent M, K, and C matrices. 

As the first step toward implementing the 
procedure using a block size m, a number of rele- 
vant matrices are defined as 

Gi = n x m complex matrix with 

orthonormal columns, that is, 

_T 

GiGi = I 

E-j = m x m Hermitian matrix 

Fi = m x m upper Hessenberg complex matrix 

Wi , X-j = n x m complex matrices 

Ti = im x im block tri-diagonal Hermitian 
matrix with blocks of size m x m 

Furthermore, a unitary matrix J p of order (p x p ) 
is next utilized to relate complex matrices 
occurring in the Lanczos method to corresponding 
matrices that are real in nature. Thus denoting S 
as a matrix of columns that are eigenvectors of H 
occurring in complex conjugate pairs, the matrix 
- -T 

S = SJ n is a real, orthogonal matrix of order n. 
Such a procedure may then be used to recast the 
block Lanczos algorithm in terms of real numbers 
effecting considerable saving in solution time. 

It may be noted in this connection that for an 
even positive integer p, the matrix J p is formed 
as p/2 replications of O 2 on diagonal, which 
in turn is defined as 



Numerical Scheme 


A = LaLa 


(5) 


in which 


*-A 



0 

k_ 


(5a) 


and Lm, are the lower triangular forms of mat- 
rices M and K, respectively. Appropriate trans- 
formation of Eq. (4) may then be effected by 
utilizing Eq. (5), yielding 


Let §1 be an arbitrary, real n x m matrix 
with orthonormal columns. Then for i = 1, 2, ..., 
the computational procedure is developed by the 
following steps. 

Step 1 . Perform matrix operation 

= HGi (for i = 1) 

A „T 

= HG-j + G-| _ 1 F i _ 1 (for i > 1) 
in which the matrix substitutions are 




(H - Yi)y = 0 

(6) 





H = i*H 

H, n x n real skew- 

in which y - 1/A, 

u) - l/i*y, and the 

matrix H is 


symmetric matrix 

expressed as 










Gi = Gi J m 

Gi , n x m real matri x 



T -T “ 


with orthonormal 


0 

-i*L M L K 



columns; that is, 

H = 

i * 

-1 -1 -T 

l K 4 i 1 * l k cl k _ 

(6a) 


Gi Gi = I 


2 


original p?-~~ 

Of POOR qimlVv 



F i = Fj, m x m real upper 

tri angul ar matri x 

Step 2 , Compute the reduced order matrix 
Ei = GiWi 

where the following substitutions are made 
-T. 


■ i * JmE i Jm 

Ej, m x m real 

skew- 

symmetric 

matri x 

= 

W-j , n x m real 

matrix 


Step 3 . Produce the matrix 

Xi = Wi - G-jE-j 

the substitution being 

Xi = i *X i J m Xi , n x m real matrix 

Step 4 . Use a standard procedure such as the 
Givens, Householder, or Gramm Schmidt method to 
obtain the "QR" factors of Xi namely Gi+i and Fi 

- - -J * 

satisfying Xi = Gi+lFi and Gi+iGi+i = I. 

Step 5 , Form the block tri-diayonal matrix 
Ti of order im 



a solution of which yields the eigenvalues and 
vectors of the system as the rth stage approxi- 
mation. Furthermore, it may be noted that T-j 

and i*Ti have the same eigenvalues that are real, 
occur in pairs, and have opposite signs. Also 
if the eigenvectors of i *f i occurring as complex 

conjugate pairs are denoted by, say, v and the 
corresponding vectors for Ti may then be obtained 
_T JT 

as v = J $ and z = J v, which are not mutually 
conjugate. The corresponding pair of real roots 
(p, -p) are Ritz values of H, the corresponding 
Ritz vectors being 

a = UiV = tfi$ 

3 = Ui z = Jiv 

where 

-T _T 

Ui = Ui J = [Gj, .... Gi ] J 

= [Si. •••» G-j J 


OF POOR QUAL5TY 


Since Qj is real, a and 0 are mutually conjugate 
so that 0 = a. 

Step 6 . Perform convergence tests using vec- 
tors computed in step 5 and matrix Fi obtained 
in step 4. If the analysis needs to be continued, 
then a selective orthogonal i zation of matrix Gi+i 

must be carried out so that its columns are or- 
thogonal i zed relative to some of the current Ritz 

vectors. ^ Thus denoting g as a column of Gi+i and 
expressing a Ritz vector a in terms of two real 
n vectors as 

a = 0 + i*# 

9 may be orthogonal i zed with respect to 0 and + and 
a new real vector § may be obtained as follows: 

g : = 9 _ ISIH e 
(e T 0) 

and 



from which the orthogonal i zed g is simply obtained 
as g - §J. 

All computations in the above procedure are 
performed in real arithmetic that has been imple- 
mented in the STARS 7 program and proves to be most 
efficient in the solution of vibration problems of 
complex, gyroscopic systems. 

Numerical Examples 

The newly implemented block Lanczos procedure 
employing real numbers (BL/R) is used to solve an 
extensive number of test cases. Such results 
are compared with solutions obtained from other 
existing similar techniques such as the block 
Lanczos technique using complex arithmetic (BL/C) 
and the SS/II methods. Because all such proce- 
dures have been implemented in the STARS program, 
it was used to perform analyses of a number of 
test cases presented here, employing a Digital 
Equipment Corp. VAX 11-750 computer. 

Spinning Cantilever Beam 

A spinning cantilever beam (Fig. 1), discre- 
tized by 12 line elements for the natural fre- 
quency analysis, has the following relevant 
properties 


Element length (a) 

5.0 


Moment of inertia (lyy) 

1/12 


Moment of inertia (Izz) 

1/24 


Cross-sectional length (A) 

1.0 


Young*s modulus (E) 

30 x 

106 

Element mass/unit length 

1/5 


Uniform spin rate (fly) 

0.33 

Hz 



and results of such analyses employing various 
procedures are given in Table 1. 

Spinning Cantilever Plate 

Figure 2 depicts a rectangular cantilever 
plate spinning along an arbitrary axis with a 
uniform spin rate Ar. A 10 x 15 finite-element 
mesh employing thin-shell elements is used to 
model the plate that has the following structural 
characteri sties 


X-side length {l%) 

10 

Y-side length ( £ y ) 

15 

Thickness (t) 

0.1 

Young's modulus (E) 

10 x 106 

Mass density (p) 

0.259 x 10*3 

Poisson's ratio ( v) 

0.3 

Number of degrees of freedom 

1056 


The plate was first analyzed for a spin rate 
= 0. 7d)i and subsequently for a resulting spin 
vector flp = 0.7u>i having components flx = = n Z ~ 

0.7 uq/VT. Results of such analyses by the three 
solution techniques are given in Table 2. 

A comparison of results presented in the two 
tables amply demonstrates the significant advan- 
tages of the present procedure. 

Concluding Remarks 

A block version of the Lanczos algorithm has 
been presented that exploits matrix sparsity and 
further performs all numerical computations in 
real numbers for the eigenproblem solution of gyro- 
scopic systems. While each solution step in the 
block algorithm is costlier than the conventional 
nonblock Lanczos method, 6 fewer steps are needed. 
The overal 1 saving in solution time is comparable 
to that effected by block multivector inverse 
iteration in place of the si ngl e-vector iteration 


process. Furthermore, this procedure is capable 
of effective determination of multiple roots in 
which the usual nonblock procedure is deficient.^ 
Also, although some experience may be necessary 
in choosing an optimum block size, a range between 
2 and 4 has been found to be effective. From 
the results presented in the two tables, it is 
apparent that the current procedure is consid- 
erably more efficient than other similar existing 
solution techniques. 
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Table 1 Results of free-vibration 
analysis of cantilever beam spinning 
at rate of 0.33 Hz (2.073 rad/sec) 


Mode 

Eigenvalue 

BL/R 

BL/C 

SS/II 

1 

3.345 

3.345 

3.350 

2 

3.602 

3.604 

3.603 

3 

16.403 

16.404 

16.404 

4 

22.491 

22.490 

22.491 

5 

44.486 

44.487 

44.487 

6 

62.258 

62.259 

62.258 

CPU time. 

20 

75 

73 

sec 





BL/R = block Lanczos procedure using 
real numbers 

BL/C = block Lanczos procedure using 
complex numbers 

SS/II = combined Sturm sequence and 
inverse iteration procedure 

Central processing unit time is for 
10 modes and frequencies 


Table 2 Natural frequencies of a spinning cantilever plate 




Natural frequency (rad/sec) 



«R = 

ftl * 0.70u>i 

= 

«R 

= 0.70un 


Mode 

149.50 rad/sec 


50 

X 

It 

» 

-< 

It 

OZ = 86.32 

rad/sec 


BL/R 

BL/C 

SS/II 

BL/R 

BL/C 

SS/II 

1 

526.69 

526.71 

526.71 

319.46 

319.46 

319.46 

2 

780.31 

780.30 

780.30 

522.79 

522.79 

522.79 

3 

1375.15 

1375.10 

1375.14 

933.32 

933.32 

933.32 

4 

1734.95 

1735.03 

1735.03 

1339.69 

1339.71 

1339.67 

5 

1791.21 

1791.25 

1791.25 

1586.76 

1586.83 

1586.83 

6 

2613.78 

2615.02 

2615.02 

2049.31 

2049.51 

2112.78 

CPU time. 

15 

52 

251 

14 

45 

255 

min 








BL/R = block Lanczos procedure using real numbers 

BL/C = block Lanczos procedure using complex numbers 

SS/II = combined Strum sequence and inverse iteration procedure 

Central processing unit time is for 10 modes and frequencies 


5 








1 Report No. 2. Government Accession No. 

NASA TM-88290 

3. Recipient's Catalog No. 

4 Title and Subtitle 

IMPLEMENTATION OF A BLOCK LANCZOS ALGORITHM 
FOR E IGENPROBLEM SOLUTION OF GYROSCOPIC SYSTEMS 

5. Report Date 
April 1987 

6. Performing Organization Code 

7 Author(s) 

Kajal K. Gupta, NASA Ames-Dryden, and 
Charles L. Lawson, Harvey Mudd College 

8. Performing Organization Report No. 
H-1404 

10. Work Unit No. 

RT0P 533-02-51 

9. Performing Organization Name and Address 

NASA Ames Research Center 
Dryden Flight Research Facility 
P.0. Box 273 
Edwards, CA 93523-5000 

11. Contract or Grant No. 

13. Type of Report and Period Covered 
Technical Memorandum 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 

14. Sponsoring Agency Code 

15 Supplementary Notes 


Prepared as AIAA Paper 87-0946-CP to be presented at the AIAA Dynamics Specialty Conference, 
Monterey, California, April 9, 1987. 


16. Abstract 


This paper describes the details of implementation of a 
general numerical procedure developed for the accurate and 
economical computation of natural frequencies and associated 
modes of any elastic structure rotating along an arbitrary 
axis. A block version of the Lanczos algorithm is derived 
for the solution that fully exploits associated matrix spar- 
sity and employs only real numbers in all relevant computa- 
tions. It is also capable of determining multiple roots and 
proves to be most efficient when compared to other, similar, 
existing techniques. 


17. Key Words (Suggested by Author(s) ) 


18. Distribution Statement 



Block Lanczos algorithm 
Eigenvalue problem 
Finite elements 
Spinning structures 
Structural dynamics 


Unclassified — Unlimited 

Subject category 39 

19. Security Ctassif. (of this report) 
IJncl assi fied 

20. Security Classif. (of this page) 
Unci assi fied 

21. No. of Pages 
7 

22. Price* 

A02 


1 For sale by the National Technical Information Service, Springfield , Virginia 22161 











